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ABSTRACT 

We present the results of near-infrared (2 and 3 /xm) monitoring of Sgr A*- 
IR with 1 min time samphng using the natural and laser guide star adaptive 
optics (LGS AO) system at the Keck II telescope. Sgr A*-IR was observed 
continuously for up to three hours on each of seven nights, between 2005 July 
and 2007 August. Sgr A*-IR is detected at all times and is continuously variable, 
with a median observed 2 fim flux density of 0.192 mJy, corresponding to 16.3 
magnitude at K'. These observations allow us to investigate Nyquist sampled 
periods ranging from about 2 minutes to an hour. Using Monte Carlo simulations, 
we find that the variability of Sgr A* in this data set is consistent with models 
based on correlated noise with power spectra having frequency dependent power 
law slopes between 2.0 to 3.0, consistent with those reported for AGN light curves. 
Of particular interest are periods of ~ 20 min, corresponding to a quasi-periodic 
signal claimed based upon previous near-infrared observations and interpreted 
as the orbit of a 'hot spot' at or near the last stable orbit of a spinning black 
hole. We find no significant periodicity at any time scale probed in these new 
observations for periodic signals. This study is sensitive to periodic signals with 
amplitudes greater than 20% of the maximum amplitude of the underlying red 
noise component for light curves with duration greater than ~ 2 hours at a 98% 
confidence limit. 

Subject headings: black hole physics. Galaxy: center, techniques: high angular 
resolution 
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Introduction 



The existence of a super-massive black hole with a mass of ~ 4 x IO^Mq at the center 
of the Galaxy has now been firmly established from monitoring the orbits of the stars in the 
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states when its flux increases by up to an order of magnitude over time scales of 1 to 3 
hours, Sgr A* is difficult to detect in its faintest states at X-ray wavelengths because of the 
strong diffu se background, and in the near-infrared be cause of confusion with nearby stel- 
lar sources (iBaganoff et al.ll2003l ; iHornstein et al.l 120071 ). Advances in adaptive optics (AO) 
technology have offered improved sensitivity to infrared emission from the location of Sgr 
A* against the stel lar background, such that observations in its faint states are now possible 
( iGhez et al.ll2005bl ). Hereafter, the IR-luminous source Sgr A*-IR will be referred to simply 
as Sgr A*, recognizing that it is likely to be coincident with the radio source of that name. 

At both NIR and X-ray wavelengths, a possible quasi-periodic oscill ation (QPO) signal 
w ith a ~ 20 min period has been report ed in light curves of Sgr A* fiGenzel et al.l 12003 



Aschenbach et al.ll2004j : lEckart et al.ll2006l ). Models that aim to produce QPO signals include 
both a class of models involving the Kepleria n orbits of 'hot spots' of plasma at the last stable 



orbit (IMeyer et al.l l2006l ; iTrippe et al. 



20071) as well as m odels with rotational modulations 



of instabilities in the accretion flow (IFalanga et al.l 120071 ) . Since the orbital period at the 
last stable orbit of a non-spinning black hole is 32 (Mb/i/4.2 x lO^M©) min, this putative 
periodic signal has been interpreted as evidence for a spinning black hole. The challenges for 
these claims are the relatively short time baselines of the observations (only a few times the 
claimed period), the low amplitude of the possible QPO activity, and the level of rigorous 
assessment of the statistical significance of the claimed periodicity. 

An alternative explanation for peaks in the periodograms seen in previous studies and 
interpreted as a periodic signal is that they are a sign of a frequency dependent physical 
process, commonly known as red noise (|Presslll978l ). The power spectrum of such a physical 
process will display an inverse power law dependence on frequency, which manifests as light 
curves with large amplitude variations over long time scales and small amplitude variations 
over short time scales. The power spectrum of any individual realization of a red noise light 
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curve will show statistical fluctuations around the intrinsic power law function, creating 
spurious peaks that can lead to an interpretation of periodic activity. Variability studies of 
other accreting black hole systems like AGNs and Galactic X-ray binaries have shown that 
their power spectral densities are consistent with red noise. Several phy sical models have 
been proposed to produce t he red noise light cu rves seen in AGNs (e.g. iLyubarskiil 119971 : 
Armitage fc Reynolds! l2003l : IVaughan et al.l 120031 ): one common model that results in a red 
noise spectrum is from fluctuations in the physical parameters, such as the gas densities 
and accretion rate, a t different radii of a turbulent magnetohydrodynamic accretion disk 
( iKataoka et al.l 120011 ) . While QPO signals have been unambiguously confirmed in X-ray 
binari es, no QPO signals i n AGNs have b een shown to be st atistically differ e nt than red 
noise (IBenlloch et al.l l200ll : IVaughanI |2005| ) . Recent work by iBelanger et al.l (Isubmittedl ) 
have also shown that a statistical analysis of X-ray light curves of Sgr A*, when including 
the contribution from red noise, show no indications of a QPO signal. 

High sensitivity and high angular resolution near-infrared observations of Sgr A* have 
been obtained at the Keck II telescope utilizing new improvements in adaptive optics tech- 
nologies to investigate the existence of a QPO signal as well as the timing properties of Sgr 
A*. These observations are described in Section [21 In order to test whether the variability 
of Sgr A* has the characteristics of red noise and to examine the possibility of a periodic 
signal, we have carried out a statistical analysis of the timing properties of the observed light 
curves that includes the possible contribution of red noise in the power spectrum to establish 
the significance of peaks in the periodograms. We find that the near-infrared variability of 
Sgr A* is entirely consistent with red noise, with no periodic signals detected on any night. 
Sections 13.2.11 details our analysis of the light curves. Lastly, in Section HI we discuss the 
implications of our results for Keplerian models of the Sgr A* flux variability and compare 
the timing properties of Sgr A* with those of AGNs. 



2. Observations and Data Reduction 



The Galactic center has been extensively imaged between 2005 and 2007 with the Keck 
II 10 m telescope using the natural and laser g uide star adaptive optics (NGS and LGS AO) 
( IWizinowich et al.l l2006l : Ivan Dam et al.l |2006| ) system and the NIRC2 near infrared camera 
(P.I. K. Matthews). For this study we include all nights of LGS- AO observations at K' (2 
/im) that had sampling of 1-3 minutes, a total time baseline of at least ~ 1 hour, and at least 
60 data points. We also include one night of NGS observations at V that satisfied the same 
criteria. As summarized in Table [H this resulted in a selection of 7 data sets with durations 
ranging from 80 min to 3 hours. 
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A deta i led des cription of LGS AO observations of the Galactic center are described in 
Ghez et al.l (l2005bl ): here, we only summarize the setup for our observations. The laser guide 
star was propagated at the center of our field and for low order tip-tilt corrections, we used 
the R = 13.7 mag star, USNO 0600-28577051, which is located ~ 19" from Sgr A*. Most 
of the images were obtained using the K' band-pass filter (Ao = 2.12 /xm, AA = 0.3 fim) 
and were composed of 10 coadded 2.8 sec exposures, for a total integration time of 28 sec. 
The remaining set of observations from 2005 July 28 was taken through the L' band-pass 
filter (Ao = 3.78 fim, AA = 0.7 /im). For five of the K' nights, the time interval between 
each image is about 50 seconds, with dithers every three minutes. K' images from 2006 July 
17 were sampled at 3 minute intervals but were not dithered. The L' observations had one 
minute sampling, and were also not dithered. The 3 min dithers affect the timing analysis 
by the presense of a spike at that frequency in the periodograms (this is well reproduced by 
Monte Carlo simulations of the effects of sampling). 

Photometry was performed on the individual imag es using the point spread function 
(PSF) fittin g program StarFinder (jPiolaiti et al.l 120001 ). The program was enhanced as 
described in iHornstein et al.l (120071 ) to include the a priori knowledge of the location of the 
near-infrared position of Sgr A* and nearby sources in order to facilitate the detection of Sgr 
A* at faint flux levels in these short exposures. To do this, for each night of observation, the 
position of Sgr A* and nearby sources was determined in a nightly-averaged image produced 
by a weighted average of individual images from that night (see Figured]). We then used the 
knowledge of the location of all the sources as fixed inputs into StarFinder to more accurately 
fit for the fiux contribution of sources near Sgr A* in the individual short exposure images. 
To compensate for seeing changes through the night, a different PSF was constructed for 
each image. We also include only images with Strehl ratios greater than 20% to minimize 
large errors in the photometry from bad seeing conditions, which resulted in dropping only 
about 10 data points out of all nights. On average, the Strehl ratio at K' was 32%, with 
the full width of the core at half-maximum intensity (FWHM) ~ 60 mas as measured from 
the relatively isolated star IRS 33N. We are able to detect Sgr A* at all times, even at its 
faintest flux levels. The gaps in the data are from technical disruptions in the observations. 

Photometr ic cal ibrations were performed relative to the list of non- variable sources from 
J2OO7I) at K' and IRS 16C (L' = 8.14 mag) and IRS 16NW {L' = 8.43 mag) 



Rafelski et al 



at L' (IBlum et al.lll996l ). The photometric error at each flux density level seen in Sgr A* 
was estimated by fltting a power law to the rms uncertainty in the flux for all non-variable 
stars in the same range of brightnesses observed for Sgr A* within 0."5 of the black hole (see 
Figure [2]). We flnd the typical dependence of the photometric error, a, on flux density, F, 
to be: a ~ 0.2F^'^ mJy. The flux measurement uncertainties are comparable for all nights 
except 2006 June 20, the night with the worst seeing. Within the range of observed Sgr A* 
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fluxes, we are on average able to achieve between 3 to 15% relative photometric precision 
for each 28 sec K' exposure. A source of systematic error in the flux measurements is the 
proximity of Sgr A* to unresolved sources, which could contribute flux. This contribution 



is only likely to have an impact when Sgr A* is faint (iHornstein et al.l 120071 ). but for the 
purpose of this variability study, this effect is likely only a systematic offset in the mean 
flux density and a source of white noise (for more details, see Section [3?T]) . For comparison 
to Sgr A*, we also analyze the light curves of the nearby stars SO-17 and SO-38 at K' and 
SO-2 at L' . SO-17 {K' = 15.5 mag) was chosen because it is spatially closest to Sgr A*, with 
a projected distance from Sgr A* of ~ 56 mas in 2006 May to ~ 48 mas in 2007 August; 
monitoring SO-17 is helpful to ensure that the variations in flux seen in Sgr A* are not a 
systematic effect of seeing or bias from nearby sources. The star SO-38 (~ 0.11 mJy, K' ~ 17 
mag), ~ 0".2 from Sgr A*, was chosen as a stellar reference because it has a similar flux 
to the faintest observed emission and given its proximity to Sgr A*, its photometry will be 
affect similarly from the unresolved stellar background. Figure [1] shows an image of this 
region and the location of the comparison sources with respect to Sgr A 



* 



In order to characterize possible effects on the photometry of Sgr A* by SO-17, we also 
use the photometry of two stars with separations and flux ratios similar to that of SO-17 
and Sgr A* when Sgr A* is faint. The two stars, S2-42 (K'=15.5) and S2-133 (K'=16.7), are 
located about 2" from Sgr A* and are separated by ~ 50 mas (Figured]). The rms variability 
of S2-42 and Sl-133 is about 5% and 15%, respectively, similar to the photometric precision 
we would have predicted based upon our power law fits to the rms stability of stars near Sgr 
A*. Thus, we can be confident that the photometry of Sgr A* at its faintest is similar to 
stars of that magnitude despite the proximity of SO-17. 



3. Results and Analysis 

3.1. Flux Distribution 

In order to characterize the range of fluxes observed from Sgr A*, we have constructed 
histograms of fluxes from each night as well as the combined histogram from all nights (Figure 
E]). Unless otherwise stated, the fluxes in this paper are observed fluxes and not corrected for 
extinctio n to Sgr A*. Where indicated, de-reddened fluxes have be en calculated by assumin g 



= 30 flMoneti et al.ll200lf ) and extinction law Ak' = 0.1 108 A flRieke fc Lebofskvl 119851 1. 
The comparison sources do not appear to be variable on the time scales probed in this study 
and have fluctuations consistent with Poisson noise. The flux distribution for the star SO-17 
is consistent with a Gaussian centered at 0.37 mJy at K', with a standard deviation of 0.02 
mJy; this suggests that we are able to reproduce the flux of SO-17 at the 5% level between 
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different nights. A bias in the photometry from the proximity of SO-17 to the hght curve 
from Sgr A* should manifest itself as a difference in the mean flux of SO-17 between 2006 
and 2007 because SO-17 moved closer to Sgr A* in the plane of the sky between the two 
years. The fact that we observe the same mean flux from SO-17 between different years is 
also confirmation that there is little bias in the photometry of either SO-17 or Sgr A*. The 
light curves of SO-17 are also stable on each night, independent of the flux of Sgr A* and 
shows no greater variance than other stars of the same brightness in the region (within 0."5 
of Sgr A*), showing that PSF fitting from StarFinder is able to able to properly account for 
the flux of both sources. 

The cumulative distribution function for Sgr A* shows that it is brighter than SO-17 
about 15% of the time at K' . The median flux density of Sgr A* is at 0.192 mJy [K' = 16.3 
mag), or a de-reddened flux of 4.10 mJy. The flux histogram for Sgr A* is not well fitted by a 
Gaussian because it has a long tail in the distribution of flux densities at high flux densities. 
However, if the tail of the distribution is excluded, the flux distribution below 0.3 mJy, is 
well fit by a Gaussian with a mean of 0.158 mJy and a standard deviation of 0.05 mJy. The 
latter is larger than nearby sources with comparable flux densities, indicating that Sgr A* is 
intrinsically variable; for example, the flux distribution for S2-133 has a FWHM of ~ 0.02 
mJy (see Figure [5]). 

The flux density distribution of Sgr A* for the L' night shows a much larger width than 
the comparison star SO-2 (Figure Sj). The distribution appears to be symmetric about the 
mean and the CDF is consistent with a Gaussian with mean flux density of 7.69 mJy with 
a = 2.81 mJy, compared to SO-2 with a mean flux density of 7.67 mJy with a = 0.37 mJy. 

While the bright flux levels (> 0.3 mJy) can be unambiguously attributed to the black 
hole, the source of emission from the location of Sgr A* when it is faint is less certain. The 
region immediately around Sgr A* is also the location of peak stellar density, which raises 
the possibility of flux contamination from either an unresolved population of stars or a very 
faint star in a close orbit at that location. We find that the faintest observed flux density 
at K' is 0.082 ± 0.017 mJy {K' = 17.2) or 1.75±0.36 mJy de- redd ened from the location of 



Sgr A* (consistent with limits observed by lHornstein et al.l (|2002| )). Comparison of the flux 



distribution of S2-133 with that of Sgr A* below 0.3 mJy shows that, for four of the six K' 
nights, Sgr A* has a larger variance than S2-133 (see figure[5]), indicating that Sgr A* is more 
variable than expected for a stellar source even at the faintest levels. On the remaining two 
nights, the variance of Sgr A* is similar to that of S2-133. On one of these two nights, 2007 
August 12, Sgr A* was fainter than 0.22 mJy for the entire duration of our observation, which 
makes this night ideal for timing analysis of Sgr A* at its faintest flux density levels. Though 
the flux distribution looks similar to a star, the structure function and the periodogram shows 
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a slightly steeper slope than expected for Gaussian noise and as compared to the stellar stellar 
comparison sources (see Sectio ns 13.2.11 and I3.2.2p. Fu rthermore, the K' — L' color for 2006 



July 17, previously reported in iHornstein et al.l (120071 ). is constant and significantly redder, 
even at its faintest on that night, than from a stellar source. At the faintest fluxes, between 
0.10 and 0.15 mJy, the mean K' — L' spectral slope, corrected for extinction, of Sgr A* 
has an average power law exponent of -0.17±0.32, compared to a slope of —0.6 ± 0.2 from 



Hornstein et al.l (120071 ). We estimate that a stellar source, which would have de-reddened 



spectral slope of 2, could contribute a maximum of about 35% of the flux to account for the 
difference in the spectral slope. This leads us to conclude that, even when the emission is 
faint, a large fraction of the flux arising from the location of Sgr A* is likely non-stellar and 
can be attributed to physical processes associated with the black hole. 



3.2. Light Curves and Timing Analysis 

Figures ini and [7| shows the resulting light curves for Sgr A* and a non-variable comparison 
source for each night of observation. While comparison sources show no significant time 
variable emission, Sgr A* shows variations on time scales ranging from minutes to hours, 
with peak emission that can be 10 times higher than during its faintest states. The emission 
peaks are time symmetric, with similar rise and fall times. 

In order to characterize the variability of Sgr A*, we have carried out the following three 
different approaches to timing analysis: (1) periodograms (2) structure functions and (3) 
auto-correlations. The periodogram analysis, presented in section [3.2. 11 is effective at pulling 
out periodic structure in light curves and therefore is optimal for assessing the presence of any 
periodicities, such as the proposed ~ 20 min QPO reported in previous experiments (e.g.. 



Genzel et al.l l2003l : lEckart et al.ll2006l ). Both the periodogram and the structure function 
can also be used to measure the underlying power spectral density (PSD), which can then 
be used to explore similarities to the variability observed in AGNs. We also compute the 
auto-correlation for each night to look for possible differences in the variability at each time 
scale between each night. 



3.2.1. Periodogram 

It is important to consider all possible sources of noise when testing for periodicity 
in light curves. While peaks in the periodograms are a good place to start searching for 
periodicity, the peaks must have significantly more power than those produced by non- 
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periodic processes to be unambiguously attributed to a true periodicity in any variable source. 
White (Gaussian) noise processes are unlikely to lead to large peaks in the periodograms 
because they contribute equal power at all frequencies. However, time-correlated physical 
processes can result in variability that is frequency dependent. One common variability 
characteristic - often seen in AGN light curves - is red noise, which can lead to spurious 
signals in a power spectrum or periodogram from a data set having a time baseline only 
a few times longer than that of the putative period, since it will show large amplitude 
fluctuations at low frequencies and small amplitudes at high frequencies. This can lead to 
relatively large stochastic peaks in the power spectrum at low frequencies, far above what 
would be expected from white noise. We emphasize that, although the term for this type of 
power law dependence of the flux variability is 'red noise', this variability arises from physical 
processes from the source and is not a result of measurement uncertainties such as Poisson 
noise, which behaves like white noise in its power spectrum. 

One of the goals in this timing analysis is to test whether a purely red noise model 
can explain the variability of Sgr A*. The PSD of a red noise light curve is a power law, 
with greater power at lower frequencies: P{f) = /~", where / is the frequency and a is 
the p ower law ind ex. For example, a = for white noise and a = 1 for classical flicker 



noise (|Presslll978f). All red no i se sin iulations in this paper were produced by an algorithm 



detailed in iTimmer &: Koenigl (Il995l ). which randomizes both phase and amplitude of an 
underlying power law spectrum and then inverse Fourier transforms it into the time domain 
to create light curves. Our procedure for producing simulated light curves is as follows: 
(1) a light curve is produced from a PSD with a specific power law slope evenly sampled 
at half the shortest observed time sampling interval, with a duration at least 10 times as 
long as the observed light curve (rounded up to the nearest power of 2 for computational 
ef ficiency of the fast F ourier transform). This length was chosen based upon the suggestion 



by lUttley et al.l (120021 ) to avoid a 'red noise leak' where power is distributed from frequencies 
lower than that sampled by the observation into observed frequencies. We find aliasing to 
be a negligible effect, because the light curves are generated at higher temporal resolution 
than the observations. ( 2) this light curve is then split into 10 non-overlapping segments to 



reduce simulation time (lUttley et al.ll2002l ). Each segment is then re-sampled at the exact 
sampling times used during the specific night that we are simulating. (3) since the simulation 
has an arbitrary fiux scale, we scale the light curves to have the same mean fiux level and 
standard deviation as that night. We also include the effects of measurement noise in the 
simulations by adding Gaussianly distributed noise to each simulated data point. We use 
our measurements of the photometric error as a function of flux densities (Section [2]) of 
non-variable stars within 0".5 of Sgr A* to account for the flux density dependence in the 
noise for each simulated data point. 
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Instead of computing the PSD, which is often used for evenly sampled data, we searched 
for periodicity by computin g a related function fo r unevenly sampled data: the normalized 
Lomb-Scargle periodogram (jPress fc Rybickilll989l ): given a set of data values hi, i = 1, . . . , N 
at times tj the periodogram is defined as: 



T 



2o"^ 1 COS "^uitj - r) ^ . sin '^uj{tj - r) 



where u is the angular search frequency, h and cr^ are the mean and variance of the data 
respectively. The constant r is an offset introduced to keep the periodogram phase invariant: 

V ■ sin 2ujt^ 

tan {2ur) = ^ — ^ (2) 

2_,j cos zutj 

Since the periodogram is normalized by the variance of the flux, a light curve consisting of 
only white noise, or equivalently, red noise with a power law a = 0, will have an average 
power of 1 at all frequencies. 

The normalized Lomb-Scargle periodogram was computed for each light curve, over- 
sampled by a factor of 4 times the independent Fourier intervals in order to increase the 
sensitivity to periods between the Fourier frequencies (Figures [7] and [H]). Assuming that the 
physical source of the variability is stationary, we averaged together the periodograms for 
the five K' nights which have durations longer than 80 minutes (Figure [9]). The combined 
periodogram excludes the 2007 August 12 night because it is less than an hour long, lead- 
ing to poor sampling at low frequencies compared to the other nights. We combined the 
periodograms by averaging the Lomb-Scargle power at linearly space frequency bins. The 
combined periodogram is consistent with red noise, except for the peak corresponding to 
the time scale of the three minute dithers. To characterize the underlying spectrum, we 
have performed Monte Carlo simulations combining red noise light curves with the same 
sampling as the data set. We tested several different underlying PSD and found that the 
combined periodogram is consistent with power law indices between 2.0 and 3.0, with no 
periodic components. This model is able to reproduce the slope of the periodogram, the 
increase in power at three minutes from dithering, and the flattening of the periodogram at 
very low frequencies caused by poor sampling at those frequencies. The 3 min peak in the 
periodogram is repoduced very well by the simulations, showing that the simulations are 
correctly accounting for the effects of sampling. Figure [9] shows the results of Monte Carlo 
simulations with power law indices 1.5, 2.0, 2.5, and 3.0. The simulations shows that the 
resulting periodograms tend to be flatter than the intrinsic PSD because the limited time 
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sampling at low frequencies results in poor sensitivity to long time scale variations character- 
istic of steeper power laws. Because the periodogram sampling is poor at frequencies below 
40 min we will ignore those frequencies in all subsequent analysis. We find the resulting 
periodograms show less variation for intrinsic PSD a > 2, which suggests that we have a 
better constraint on the lower limit than on the upper limit of our estimate for the slope 
of the PSD. While red noise models with values of a between 2.0 to 3.0 appear to be con- 
sistent with the average periodogram, we will use a = 2.5 for the simulations in this work 
to provide a baseline for comparison. More light curves will be necessary to determine a 
reliable intrinsic PSD of Sgr A* (if it does not vary between nights). Where appropriate, we 
have also run simulations with a range of a values to investigate its effects on the statistical 
significance. 

Although there is no evidence for QPO activity in the combined periodogram, we also 
test the case of a transient QPO phenomenon in each night. We therefore used our best fit 
of a = 2.5 for the power law of the combined periodogram to test each night for statistically 
significant deviations from a purely red noise model. Our criterion for a statistically signif- 
icant peak in the periodogram is that its power must be above the 99.7% (Sex) confidence 
interval from a Monte Carlo simulation with 10^ realizations of red noise light curves. This 
method of e stablishin g the significance of peaks in the periodogram is similar to the one 



proposed by IVaughanl (120051 ) for evenly sampled data, but modified here to account for our 
unevenly sampled data set by using the Lomb-Scargle periodogram instead of the Fourier 
transform. The resulting periodograms are shown in Figure [Hi The individual periodograms 
show peaks at low frequencies, but these appear to be consistent with red noise, with no 
peaks having power greater than the 3a threshold derived from the Monte Carlo simulations. 
No significant peaks were found when we repeated the same procedure for PSD power law 
indices between 1.0 and 3.0. 

The periodograms of the comparison sources are much flatter than that of Sgr A*, as 
expected for white noise sources. However, there is some power at very low frequencies (< 
40 min), which may be the result of some small time correlated systematic in the observation 
such as variations in seeing or AO performance over a night. Light curves and periodograms 
from 2006-06-21 of the three comparison sources SO-17, SO-38, and S2-133 are shown in 
figure Uni The periodograms of these sources from the other nights of observations are also 
comparably flat. 

By adding an artificial sinusoidal periodic signal to the red noise simulations, we can 
address our sensitivity for detecting QPOs in the presence of red noise. In order to test this, 
a periodic sinusoidal signal was added to a simulated light curve with a red noise PSD slope 
of 2.5, with an amplitude that is a fraction of the maximum signal in the underlying red 
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noise light curve. This fractional amplitude of the periodic signal was increased until the 
periodogram showed a peak at that frequency above the previously determined 3cr threshold 
over 98% of the time. For example, for the time sampling and duration corresponding to 
the observations on 2006 May 3, we find that we are able to detect a 20 min periodic signal 
with an amplitude that is at least 20% of the maximum flux density (Figure [TTj) . The 
sensitivity for the detection of a QPO increases with frequency because the underlying red 
noise component has less power at higher frequencies (e.g. a 10 (40)-minute QPO signal will 
result in power greater than expected from red noise when its amplitude is 5% (30%) of the 
maximum flux density). Our sensitivity for the detection of a periodic signal is similar for 
the other observed light curves with similar durations (> 2 hours). 

We note that attempting to remove the broad, 40-100 minute emission peaks with a low- 
order functional fit to increase sensitivity to periodic signals will introduce a statistical bias, 
because that will only remove some combination of low frequency power without actually 
removing the entire red noise component. Because the light curves are consistent with 
a red noise process at all time scales, including the longer broad maxima, any statistical 
analysis must be performed without fi rst modifying the light curves. For example, the light 
curve reported by iGenzel et al.l (120031 ) . the first used to argue for the possible presence of a 
periodic variation, shows no significant periodicity when subjected to the red noise analysis 
de scribed here. Figure [T2] shows the analysis of this night using the light curve re-extracted 
by iMeyer et al.l (120061 ) . A sig nificant signa l was found when a low order polynomial fit to 
the light curve was removed in lMeyer et al.l (120061 ). However, this significance may be biased 
by the removal of a number of low frequency components from the light curve. 



3.2.2. Structure Function 



The first order structure function is often used to determine th e time scale and the 
intrinsic vari ability in AGN light curves (e.g., ISimonetti et al.l Il985l : iHughes et al.l Il992 



Paltanil Il999l ). Here, we use the structure function to determine the PSD slope for each 
night because it is more straightforward to calculate errors in the structure function. For a 
set of measurements, s{t) at times t, the first order structure function V{t) is defined as: 



V{t) =< [s(t + r) -s(t)]' > 



(3) 



Because this data set is unevenly sampled, we determined the structure function by cal- 
culating [s{t + r) — s(t)]^ for all possible pairs of time lags and place them into bins with 
widths of one fourth of the lag times. Variable bin sizes were chosen in order to more evenly 
distribute the number of points in each bin, since there are many more samples at small lags 
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than at large lags. The median lag time in each bin is assigned to be the lag time for that 
bin, while the average of the V{t) values in each bin is used as the value of the structure 
function at that lag. The error associated with each bin is ounl \J ^Un-, where oun is the 
standard deviation of the V^(t) values in the bin and Nun is the number of points. We do 
not consider lags with fewer than 5 points in that bin. The resulting structure functions are 
presented in Figure [T3l 

Structure functions typically show a power law portion and two plateaus - one at the 
time scale of the noise and one plateau at lags longer than the variability time scale of 
the underlying physical process. The power law portion of the structure function is a mea- 
surement of the variability of the process. If the underlying process is stationary, then the 
logarithmic slope, /3 (where Vij) oc of this region can be related to the power index, 
a of the PSD iP oc /^"). Because the light curves in this study have finite sampling and 
duration, this relationship must be determined separately for each night. To convert the 
measured structure function slope into the PSD slope, we measured the slope of an average 
of 10^ structure functions (/?) from simulated light curves for a number of PSD slope values 
[a). The relationship between the measured structure function slope and the PSD slope is 
linear for values of a between 1.0 and 3.5 (Figure fT6i) . These linear fits are used below to 
convert between the two slopes. We will use the structure function to compute the PSD 
because it is more straightforward than fitting the periodogram, given the limitations of the 
time sampling. 

The structure function for Sgr A* shows a power law portion at periods greater than 
about 1 minute and a flattening at periods greater than about 40 minutes. This flattening is 
probably not a sign of a turnover in the intrinsic variability time scale, but is rather from the 
poor sampling at longer periods, as shown by large variations in the value of the structure 
function from Monte Carlo simulations of red noise using the same time sampling. This 
suggests that Sgr A* is likely to be variable at all the time scales probed in this study. Fits 
were made to the power law region of the structure function, typically up to about half or 
one-third of the total time of observations for each night, and are summarized in Table [2l 
The best fit power law index for the nights of 2006 May 3 and 2006 June 20, and 21 has 
/? = 1.4 ± 0.1 [a = 2.7 ± 0.1), between 1 to 20 minutes, consistent with the fit of a ~ 2.5 
for the slope of the combined periodogram. The other three nights have slopes ranging from 
P ~ 0.8 to /9 ~ 0.3. The night with the flattest slope (2007 August 12) is also the night where 
the emission from Sgr A* was faint for the duration of the observation. Figure shows that 
though the slope, j3 = 0.26 ± 0.04, is flatter than on other nights, it is significantly steeper 
than that of SO-17 {(3 = 0.04±0.04), S2-133 {(3 = 0.11 ±0.05), and simulations of white noise 
(/3 = —0.01), which shows that the emission is not from a constant background source. The 
differences in slope between different nights may be attributable to the increase in photon 
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noise at these faint fluxes (see Section Hj). For all nights, the SO-17 structure function is 
also flat, with a slight rise at longer time scales (also from poor sampling). The structure 
function of Sgr A* for the L' night on 2005 July 28 has a best fit slope /3 = 0.77 ± 0.18, 
corresponding to a = 1.84 ±0.27, which is within the range of observed K' power law slopes 
(Figure [T5|) . 



3.2.3. Auto -correlation 

To compare the Sgr A* variability between each night, we computed the auto-correlation 
function for each observation. The auto-correlation is defined as: 

AC{t) =< [sit + r)s(t)] > (4) 

where s{t) is the measured flux at time, t, and r is the time lag. Because the light curves 
are unevenly sampled, we computed the auto-correlation at all possible lags and then placed 
the pairs in bins, similar to our method for the computing of the structure function. We 
find that the correlation coefficients at time scales up to 40 minutes are consistent across 
all nights. At time lags greater than 40 minutes, the auto-correlations appear to diverge 
from each other. This is mostly likely from poor sampling at those frequencies. From Monte 
Carlo simulations of red noise with the same time sampling as the data set, we find that the 
lags greater than about 40 minutes have large variations in auto-correlations coefficients; we 
therefore defer interpretations of lags greater than 40 minutes until more data is available 
to better sample these time scales. 



4. Discussion 

Our analysis shows that the Sgr A* near-infrared light curves observed in this study are 
entirely consistent with red noise with no significance evidence for a QPO signal near 20 min 
(or any other time scale). In addition, power is seen at the shortest time scales measured 
in the periodogram (2 min), well below the orbital period at the last stable circular orbit 
of a spinning black hole (~ 5-30 min). This indicates that the high frequency variability is 



Mever et al. 


2006; 


Trippe et al. 


2007) 



originate from small areas of the accretion flow. In addition to power at high frequencies, 
the light curve for 2007 May 18 shows a strong spike in flux at about 10 min with a duration 
of ~ 5 min, and a rise time of only ~ 3 min, which also suggests that regions of the accretion 
flow responsible for the variability can be as small as 0.4 AU. 
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It is important to note that the significance threshold used in this study tests whether 
the periodogram power of peaks exceeds 3a at a given frequency. It does not take into account 
that we are scanning over a large range of frequencies. If one searches enough frequencies, 
a peak of arbitrary power will appear at a random frequency. If a peak in the periodogram 
was found in this study that may indicate a periodic signal, then additional simulations 
should be performed to determine its statistical significance by accounti ng for the range of 
freque ncies probed. This is referred to as the false alarm probability in iHorne fc Baliunas 
(119861 ) and has been calculated analytically for the case of white noise. In the case of red 
noise, the additional simulations would also need to account for the fact that larger peaks 
are more likely at low frequencies (IMeyer et al.ll2008l ). 

The timing analysis of Sgr A* in the infrared may be a way to compare the signatures 
of physical processes operating in the accretion flow of Sgr A* to those of AGNs. Because 
AGN variability studies are usually performed on X-ray observations, the ideal comparison 
of the variability of Sgr A* to that of AGNs should be done at X-ray energies. However, 
X-ray emission from the inner accretion flow of Sgr A* can only be currently detected above 
the diffuse background during large flaring events. Simultaneous multi-wavelength studies of 
Sgr A* have shown that every d etectable X-ray flar e shows corresponding pe aks in the near- 
IR, with no apparent time lag (lYusef-Zadeh et al.l l2006l : lEckart et al.l 120061). The spectral 



slope is also similar at X-ray and infrared wa velengths {Si, oc u : iBelanger et al.l 12005 



Yusef-Zadeh et al.ll2006l : iHornstein et al.l 120071 ) . which supports a model in which the X-ray 



photons are the result of synchrotron self Compton scattering of lower energy photons from 
the population of electrons that produce synchrotron emission seen in the infrared. If the 
infrared and the X-ray emissions are produced by the same population of electrons, then 
our infrared measurements can be used as a proxy for the X-ray variability; we flnd that 
the periodograms of the infrared light curves of Sgr A* are in fact consistent with having 
intrinsic PSD slopes sirnilar to those reported for AGNs at X-ray waveleng ths {a ~ 1—3) (e.g. 



Lawrence fc Papadakislll993l : iKataoka et al.ll200ll : iMarkowitz et al.ll2003l ). Our observations 
are consistent with Sgr A* being variable at all time scales between ~ 2 minutes to ~ 
1 hour with no changes in PSD slope between these time scales during the observations 
period on any single night. This steeply rising power law cannot be extended to arbitrarily 
low frequencies because the integrated power would diverge. For example, AGNs typically 
show either one or two breaks in their power spectrum, with rece nt evidence suggestin g 
that these breaks scale with the mass and accretion rate of the AGN (IMcHardy et al.ll2006l ). 
The lack of a break in the Sgr A* periodogram is not unusual because the br eaks seen in 
AGNs are at time scales on the order of days or weeks (IMarkowitz et al.ll2003l ). Using the 
relationship between t he break frequency, the black hole mass, and accretion rate found by 



McHardy et al.l (120061 ) for stellar mass black holes and AGNs, the predict break frequency 
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for Sgr A* is ~ 1 X 10^ days. However, a direct comparison to AGNs is difficult at this time 
because the physical scales involved are so different; the cores of AGNs in which variability 
is usually measured have sizes on the order of parsecs, whereas Sgr A* is unresolved in this 
study at a resolution of 65 mas, or ~ 540 AU. The accretion luminosities of AGNs with well 
measured timing properties are substantial fractions of their Eddington luminosity compared 
to the 10~^ Eddington luminosity for Sgr A*, which may make the comparison inapplicable 
as well. Nevertheless, the variability properties of Sgr A* are valuable additions to the 
broadband spectrum since any models that can produce the observed spectrum should also 
be able to account for the timing characteristics. 

While it may appear from the results in Section 13.2.21 that the slope of the power spec- 
trum is changing between nights, the red noise simulations show that, as the amplitude 
of the red noise signal decreases, photon noise will increasingly dominate the measurement 
of the slope, thus biasing the measurement of the true slope of the underlying red noise 
process. To determine when photon noise becomes important to the slope measurements, 
we performed Monte Carlo simulations of light curves with time sampling from 2008 Au- 
gust 12 and assuming an intrinsic a = 2.5 for various values of the ratio of the standard 
deviation of the red noise signal to the standard deviation of the expected measurement 
noise (cr^ed noise/ o'whue noise)- Average values of a were extracted from 200 simulated light 
curves at each value of ared noise/ CTwhite noise from 0.1 - 100.0 (Figure [T71). We find that as 
o'red noise/ o'white noise bccomcs Smaller, the measurcd value of a becomes flatter. The simu- 
lation showed that a measurement of a = 1.0 can arise even from a red noise process with 
a = 2.5 if ared noise/ o'white noise = l-l- Wc find that the measured value of a will asymp- 
tote to 2.5 when (Xred noise/ o'white noise > 7.9. This rcsult suggests that the measurement of 
a = 1.06 ±0.08 for 2007 August 12 can be explained as having the same intrinsic PSD slope, 
a, as the other nights, but made flatter by a more prevalent white noise component. Note 
that this is an effect of a white noise source on top of a red noise signal, which can be caused 
simply by the measurement noise at faint flux densities or the presence of an unresolved 
stellar source. Figure [17] shows the relationship between the measured value of a and the 
ratio of standard deviation of the observed light curves to the expected noise level for the 
corresponding mean flux density {o signal /onoise) for all K' nights. The measured value of a 
appears to asymptote on nights when the red noise amplitude is stronger. This relationship 
is consistent with that found in the Monte Carlo simulations, where the nights with lower 
values of o signal / o noise havc flatter red noise spectra than those on the nights with greater 
fractional signal variance. These results show that the differences between the measured 
slopes for different nights may be explained by the effect of measurement noise on light 
curves with faint fluxes. 

We also considered how a slightly time correlated source would affect the power law 
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slope measurements by including a time correlated source with a small a into the simulations 
and varying its signal strength as the above treatment with white noise. This is motivated 
by the fact that there is a slight slope to the structure function of our comparison source 
S2-133. For example on 12 August 2007, S2-133 has a structure function power law slope 
of f3 = 0.11 ± 0.05 corresponding to a = 0.7 ± 0.2. We find that the addition of such a 
background source increases the measured Sgr A* PSD power law slope a by a maximum of 
0.1, which is about the same level as our fitting error. The effect of a shghtly time correlated 
background source thus appears to be a source of systematic error, but one which does not 
significantly affect our results. 

A notable feature of the light curves in these observations is that there appears to 
be a lower limit, or 'floor' to the light curves from each night. This lower limit is not an 
observational constraint, because we are able to rehably detect stars within 0".5 of Sgr A* 
to K' Ri 19 mag compared to the faintest Sgr A* observation at K' 17 mag. A question 
to be investigated in future studies is whether this represents a 'quiescent' state of Sgr A*, 
characterized by lower variability at that flux level or is a faint star confused with Sgr A*. If 
the low flux density state is a physically distinct state, then the large flux variations would 
be the result of additional emission of energy above this quasi-steady state. To examine the 
nature of this faint state, we can look to the peak of the distribution of flux densities fainter 
than ~ 0.3 mJy for each night. We find that the fiux density of the peak of this distribution 
shifts slightly between different nights of observations. However, this shift appears only to be 
marginally significant given the lower photometric accuracy at the fainter fiux levels. Before 
physical models can be applied to the faint state, a more detailed analysis of the light curves 
is necessary to establish whether this state exists. If the faint state of Sgr A* can be modeled 
as arising from the steady accretion of mass, then the slight drift of this mean flux level may 
be explained by a gradual change in the accretion rate. Another explanation for the drift in 
the mean Sgr A* flux density is that it is due to very long timescale flux variation from the 
red noise power law behavior of the flux as in AGNs. 

It is clear that Sgr A* can now always be detected in the infrared using LGS AO and 
that even at low flux densities, the K' — L' color of the source appears to be redder and have 
PSD slopes steeper than a star, but the measurement errors do not rule out the possibility 
of a faint stellar component in the near-infrared fiux from the location of Sgr A*. More 
observations during both the faint and bright states of Sgr A* are necessary to confirm 
whether there truly exists a quiescent state with different timing behavior. Since Sgr A* 
is much redder than any stellar source, observations at L' would be much less affected by 
stellar contamination than at K', thus the test for a quiescent state may be more conclusive 
at L' . Note that the existence of a quiescent state does not affect the timing analysis for 
periodic signals because the mean fiux density is subtracted from a light curve to remove 
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the zero frequency power before creating the periodogram. 

While no QPO signal was detected in this study, this does not preclude the possibility 
for a transient periodic signal in other observations. In addition, a periodic signal with an 
amplitude below our detection threshold would also be difficult to detect unless the signal 
is consistently at the same frequency such that it would become statistically significant in 
a combined periodogram. More observations of Sgr A* would be helpful to investigate such 
low amplitude QPO signals. 



5. Conclusion 

We have obtained 7 near-infrared light curves of Sgr A* using NGS and LGS AO and 
through Monte Carlo simulations, test for a periodic signal in the presence of red noise with 
unevenly sampled data. We find no statistically significant QPO signals in any of the Sgr A* 
light curves observed in this study; the variability appears to be consistent with a physical 
process having a power spectrum index a ~ 2.5 in the averaged periodogram, with variations 
between a = 1.0 to 2.8 between different nights. A determination of the true PSD power law 
slope of Sgr A* will require at least double the number of light curves observ ed in this study 



so that more sophisticated sp ectral estimation techniques can be used (e.g. |Priestleylll981 



Papadakis &: Lawrence! Il993l ). While more data will be necessary to reliably determine the 
intrinsic PSD slope, it is unnecessary to invoke a periodic signal to describe the variability 
of the light curves observed in this study. We have also shown the importance of preforming 
statistical tests for QPO signals using the entire Sgr A* light curve. Removing broad large 
amplitude peaks eliminates some combination of low frequency power, artificially boosting 
the significance of mid-range frequencies without actually removing the underlying red noise 
component. 

The mechanism producing the infrared variability in Sgr A* in these observations may 
be related to the processes occurring in AGNs or accreting black hole binaries, based on the 
similarity of their timing properties. Detailed model comparison between Sgr A* and other 
accreting black hole systems will require continuous observations over longer time scales and 
at faster time sampling to detect possible low and high frequency breaks, respectively, in the 
power spectrum. 
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Fig. 1. — Left: a K' image from 2006 May 3 of the central 0."5 around Sgr A* with 

a logarithmic intensity scale so that the faint sources can be more easily seen. Sgr A* 

(K' = 15.8 in this image) is in the center of the image along with the comparison stars, SO-2, 

SO-17 and SO-38. The image is oriented with north up and east to the left, with offsets in 

projected distance from Sgr A*. Right: image from the same night of the pair of comparison 

sources, S2-42 {K' = 15.5) and S2-133 {K' = 16.7), with a flux ratio similar to that of SO-17 

and Sgr A* when Sgr A* is faint. These two stars also has a similar separation in the plane 

of the sky as SO-17 and Sgr A* (~50 mas). 
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Fig. 2. — Plot of power law fits to the rms fiuxes for stars within 0."5 of Sgr A* with 
brightnesses within the range of observed brightness variations in the K' light curves from 
Sgr A*. The photometric noise properties are similar for each night. 
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Fig. 3. — Top: Sgr A* flux distribution for eacli of tlie 6 K' nights (left), the combined fluxes 
from all nights (middle), and the inverse cumulative distribution function for the combined 
fluxes (right). Middle: The corresponding plots for SO-17, the non- variable comparison 
source used in this study. The total histogram of fluxes observed from SO-17 is consistent 
with a single Gaussian. Bottom: the flux distribution for S2-133, a K' = 16.7 magnitude 
star for comparison with the faint states of Sgr A*. 



-24- 




Fig. 4.— Left - Sgr A* and SO-2 flux distributions on the night of 2005 July 28 at ll . Right 

- the corresponding cumulative distribution function for Sgr A* and SO-2. Both CDFs are 
consistent with Gaussian distributions. 
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Fig. 5. — The flux distribution of Sgr A* (solid) of fluxes less than 0.3 mJy and the reference 
source, S2-133 (dotted) for each night along with the best fit Gaussian. On four of the six 
K' nights, the width of the Sgr A* flux distribution is greater than the reference, while in 
the remaining two nights, the widths are comparable. Labels are units of mJy. 
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Fig. 6.— 



Sgr A* light curves (black) at K' with the star SO-17 (blue) for comparison. 
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Fig. 7. — Top: The light curve of Sgr A* (black) and comparison source SO-2 (grey) at 
V on 2005 July 28. Bottom: Normalized Lomb-Scargle periodograms of Sgr A* on this 
night (black). Also plotted is the 3 a significance threshold determined from Monte Carlo 
simulations of red noise with a power law index, a — 2.5 (dotted). The average of 10^ 
periodograms with a = 2.5 at the same time sampling is also shown for comparison (dashed 
hne) . 
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Fig. 8. — Normalized Lomb-Scargle periodograms of the Sgr A* light curves (black) at K'. 
Also plotted is the 3 a significance threshold determined from Monte Carlo simulations of 
red noise with a power law index, a — 2.5 (dotted line). The average of 10^ periodograms 
with a — 2.5 at the same time sampling is also shown for comparison (dashed line). 
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Fig. 9. — Top: averaged periodogram from 5 nights of K' observations (black line) with 
durations greater than 80 min (resulting in the exclusion of 2007 Aug 12). Colored lines are 
the average of 10'^ Monte Carlo simulations of combining data sets with the same sampling, 
for three different red noise power laws. The dotted blue line corresponds to the 3 a threshold 
of power for red noise with a = 2.5. The large spike at 0.0056 Hz (3 min) is an artifact 
resulting from the regular dithers in the observations. Bottom: plot of the number of points 
contributing to the averaged periodogram at each frequency. 
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Fig. 10. — Top: Light curves of the stellar comparison sources from 2006-06-21. The light 
curve of SO-38 is offset by 0.18 mJy for clarity. Bottom: The corresponding periodograms 
of the comparison sources. In sohd gray is the 3 a threshold for white noise (a = 0). 
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Fig. 11. — Top: Simulated light curve with the same time sampling as 2006 May 03 with a 
20 minute periodic signal with an amplitude 15% of the maximum flux seen that night. Error 
bars in the light curve are based upon the noise properties of that night. Bottom: The 
corresponding periodogram of the simulated light curve (solid hne) along with the average 
(dashed) and 3a threshold of red noise power from the Monte Carlo simulations (dotted 
line). We can recover the signal at 20 min (vertical dotted hne) at above our established 
confidence threshold. 
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Fig. 12. — Left: T he de-reddened light curve of Sgr A*, from 2003 June 16, as extracted by 



Meyer et al.l (120061 ) (filled circles) along with a comparison source, called SI in that paper 
(squares, offset by 6 mJy for clarity). Right: The periodogram of the Sgr A* light curve 
(solid) with the mean power (dashed) and 3a thresholds set by Monte Carlo simulations of 
red noise with a power law slope of a = 2.5 (dotted). 
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Fig. 13. — The structure functions for five of the K' nights for Sgr A* (sohd) and the 
comparison source, SO-17 (dashed). Error bars are calculated from a/^/N, where cr is the 
standard deviation and the number of points in the bin. Error bars for SO-17 are of 
comparable sizes, but are omitted for clarity. 
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Fig. 14. — The K' structure function of Sgr A* (black), S2-133 (blue), and simulations of 
white noise (green) for 2007 August 12.. 
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Fig. 15.— The structure function of Sgr A* (black), SO-2 (blue) at L' on 2005 July 28. The 
best fit power law between 3 to 40 min is also plotted. The best fit power law for Sgr A* is 
/3 = 0.77 ± 0.18, corresponding to a = 1.84 ± 0.27 (see Section [3221 for details). 




Fig. 16. — The relationship between intrinsic PSD power index and the measured power law 
of the structure function as simulated for the time sampling on each night. The line fits are 
made for PSD power law slope a > 1.0 in the linear region of the relationship between the 
intrinsic PSD power law and the structure function power law. 
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Fig. 17. — Plot of the fitted power law slope a as a function of the ratio between the 
standard deviation of the light curve (signal) and the photon noise expected at the mean 
flux density for each K' night (solid black squares). Open grey circles shows the results of 
MC simulations with the same time samphng as 2007 August 12, showing that as the red 
noise signal becomes comparable to the measurement noise, the measured slope of the power 
spectrum becomes flatter than the intrinsic value of o; = 2.5 used in the simulation. 
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Table 1. Summary of Observations of Sgr A* 



Date 


Filter 


Start Time 


End Time 


Nobs 


Median Samp. 


Duration 


Dithered'' Strehl*^ 


FHWM^'Mean Flux'^ 




S0-17phot 


''S2-133phot'* 


(UT) 




(UT) 


UT 




(sec) 


(min) 




(%) 


(mas) 


(mJy) 


(mJy) 


% 


% 


2005 July 28^= 


L' 


06:10:52 


09:09:23 


144 


69 


180 


No 


70 


81 


4.74 


1.7 






2006 May 03 


K' 


10:54:30 


13:14:12 


116 


51 


140 


Yes 


34 


60 


0.29 


0.15 


7 


13 


2006 June 20 


K' 


08:59:22 


11:17:54 


70 


51 


79 


Yes 


24 


71 


0.24 


0.18 


7 


22 


2006 June 21 


K' 


08:52:26 


11:36:53 


164 


51 


164 


Yes 


33 


61 


0.19 


0.07 


4 


17 


2006 July 17' 


K' 


06:45:29 


09:54:03 


70 


147 


189 


No 


35 


59 


0.15 


0.05 


5 


22 


2007 May 18 


K' 


11:34:10 


13:52:39 


77 


51 


84 


Yes 


36 


59 


0.29 


0.14 


7 


10 


2007 August 12 


K' 


06:67:10 


07:44:38 


62 


51 


57 


Yes 


33 


58 


0.16 


0.02 


4 


7 



^When the observations are dithered, the time interval between dithers is 3 min. 

''Average for the night. Strehl ratios and FWHM measurements made on IRS 33N for all K' data and IRS 16C for L' data 
'^Flux values are observed fluxes, not de-reddened. 

''Photometric precision for the two comparison sources SO-17 and S2-133, with mean fluxes 0.38 mJy and 0.15 mJy respectively. 

<=NGS mode <^ 
'Previously reported in iHornstein et al] ll2007h 
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Tablc 2. PSD fits from the stnictiirc function 



Date 


Fit Range 






(UT) 


(min) 






2006 May 03 


1-20 


1.37 ±0.09 


2.74 ±0.13 


2006 June 20 


1-20 


1.40 ±0.08 


2.79 ±0.12 


2006 June 21 


1-15 


1.42 ±0.18 


2.79 ±0.27 


2006 July 17 


1-40 


0.71 ±0.09 


1.80 ±0.16 


2007 May 18 


1-20 


0.79 ±0.11 


1.95 ±0.20 


2007 August 12 


1-30 


0.26 ±0.04 


1.06 ±0.08 



^Slope of the power law portion of the structure func- 
tion, V{t) oc 

^Slope of the PSD, P(/) oc 



